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Abstract 

Theoretical approaches of electric transport in correlated molecules usually consider an extended 
molecule, which includes, in addition to the molecule itself, parts of electrodes. In the case where 
electron correlations remain confined within the molecule, and the extended molecule is sufficiently 
large, the current can be expressed by means of Laudauer-type formulae. Electron correlations are 
embodied into the retarded Green's function of a sufficiently large but isolated extended molecule, 
which represents the key quantity that can be accurately determined by means of ab initio quantum 
chemical calculations. To exemplify these ideas, we present and analyze numerical results obtained 
within full CI calculations for an extended molecule described by the interacting resonant level 
model. Based on them, we argue that for narrower band (organic) electrodes the transport proper- 
ties can be reliably computed, because the extended molecule can be chosen sufficiently small to be 
tackled within accurate ab initio methods. For wider band (metallic) electrodes, larger extended 
molecules have to be considered in general, but a (semi-)quantitative description of the transport 
should still be possible in the typical cases where electron transport proceeds by off-resonant tun- 
neling. Our numerical results also demonstrate that, contrary to the usual claim, the ratio between 
the characteristic Coulomb strength and the level width due to molecule-electrode coupling is not 
the only quantity needed to assess whether electron correlation effects are strong or weak. 
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I. INTRODUCTION 



Intensive experimental and theoretical work done in the last decade led to significant 
advances in the field of electric transport through molecules and nanostructures 
From the theoretical side, ab initio approaches — mostly based on ground state density- 
functional theory (DFT) calculations combined with the Keldysh nonequilibrium Green's 
functions (NEGF) — contributed to a qualitative understanding of some trends impor- 
tant for molecular electronics 3|, \^ . However, intrinsic limitations of the DFT-approaches 
make it impossible to quantitatively describe electronic transport in correlated molecular 
systems. There is a broad consensus in attributing the rather large quantitative discrep- 
ancies between the calculated and measured currents to the notorious fact that the DFT 
substantially underestimates the excitation energies, in particular the HOMO-LUMO gap. 
In order to more appropriately account for electron correlation effects on nanotransport. 



several studies, which also used the 
the so-called GW approximation 



NfEGF, focused their attention on other methods, e. g ., 



10|, a method proposed long time ago for metals 



The theory of electric transport through correlated molecular systems (i. e., which cannot 
be described within a single-particle picture) is confronted with problems, which are not 
encountered in the absence of electron correlations. A basic issue is that of partitioning 



the real system into electrodes and the molecule [12] . In the uncorrelated case, the electric 
current can be computed exactly within the Landauer approach. Being exact, the result 
is partition- independent, and therefore the specification of the electrode- molecule interfaces 
used in calculations is not important. For physical reasons, the current through a correlated 
molecular system should also be partition-independent. Excepting for a few models, the 
calculations on the transport in correlated systems are not exact. The various approximation 
schemes employed do not necessarily satisfy this partition invariance. 

We are not aware of a systematic study addressing this issue. In most cases, the reported 
results are from calculations carried out for the largest extended molecule that can be tackled, 
wherein besides the molecule itself the largest possible electrode portions linked to it are 
included. The main concern of the most elaborate theoretical approaches of the correlated 
molecular transport developed so far is not the partition invariance or, alternatively, the 
independence on the size N of the extended molecule. One is merely interested whether the 
employed approximate approach is current conserving or not, i. e., whether the calculated 
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steady-state current does depend or not on the position in the extended molecule of the 

15|. 

To summarize, a reliable approach of the transport in correlated molecular systems 
should not only yield a steady-state current that does not depend on the transverse areas 
at . . . L^, L2, Li, M, i?2, -R3, • • ., but also on the way in which the extended molecule is 
delimited by the various rectangles, as schematically depicted in Fig. [H 



cross-sectional area through which it flows 
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FIG. 1: Schematic representation of a typical two-terminal setup for molecular transport. A 
molecule M is linked to two left L and right R electrodes. In a correct approach, the steady-state 
current is site-independent and does not depend on the way in which the total system is partitioned 
into the central part (extended molecule) and the electrodes. In the figure, the different partitions 
are depicted by the dashed-line frames. 

The remaining part of the paper is organized in the following manner. In Sect. [Tll we 
present the theoretical framework, and in Sect. IIIII the rationale for the concept of a suf- 
ficiently large extended molecule. The interacting resonant level model, which is used in 
all numerical calculations of this paper, is described in Sect. HVl Next, in Sect. |V] we re- 
port numerical results for the linear conductance. The case of narrower band electrodes is 
examined in Sect. IV At Sect. IV Bl is devoted to the case of wider band electrodes. To be 
specific, for these two cases we shall use the terms organic and metallic electrodes, respec- 
tively. In Sect. IVIt we present a simple renormalization scheme, which works surprisingly 
well and yields results that exhibit a very weak size dependence even for unrealistically 
strong Coulomb contact interactions. In Sect. IVIII we discuss the relevance of numerical 
results for realistic systems, and in Sect. IVIIII we present the conclusions. 
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II. THEORETICAL FRAMEWORK 



Within the Keldysh nonequihbrium Green's function framework, it is possible to express 



the steady-state current I through a correlated molecu 
trons are supposed to be uncorrelated, by the formula 



16 



inked to electrodes, wherein elec- 



18| 



I='^j dsTr [UlTi^ - fnTR) (G^ - G") + (r^ - r^) G<] . (1) 

Above, the obvious e-dependence of the integrand has been omitted for brevity. In Eq. ([T]), 
G^, G", and G^ stand for the nonequihbrium {V ^ 0) exact retarded, advanced, and the 
lesser Green's functions of the extended molecule, respectively, and and describe 
the coupling of the extended molecule to the left (L) and right (R) electrodes, respectively. 
The bold symbols are used to denote matrices, whose indices label the basis states of the 
extended molecule. These are site indices for the case considered below; in a realistic ab initio 
calculations, they would specify the molecular orbitals characterizing the extended molecule. 
As usual, each electrode is assumed to be in equilibrium; fi^R represent Fermi distribution 
functions characterized by chemical potentials /iL,_R, which are imbalanced (/xl — fJ.R = eV) 
due to an applied voltage V responsible for the current flow. 

The Green's functions G*", G", and G^ entering Eq. ([T]) are those of the extended 
molecule coupled to the electrodes. To com put e these nonequihbrium Green's functions one 



usually resorts to the Dyson equations 17|- 
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G> = (1 + G^S^) G> (1 + S'^G") + G''S>G" ^ G^S>G" , (3) 

which relate the exact Green's functions (without subscript) to those of a noninteracting 
(uncorrelated) reference (subscript 0) by means of the self-energies S's. Uncorrelated elec- 
trons can refer either to free electrons or described within the SCF approximation. The 
expression at the right of the arrow in Eq. ([3]) applies to the case of a steady-state flow 



considered here 



12| . The above nonequihbrium Dyson equations can be written in a more 



compact form by means of the block matrices [li 

/ yt -y< \ 

(4) 
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where the causal [t), anti-causal (t), and the greater (>) matrix elements are related to the 
earlier specified ones by 





G = G - G^ = G^ - G , 


(5) 




G = G - G^ = G^ - G = (G )' . 


(6) 


Dyson equations 


equivalent to Eqs. ([2]) and ([3]) relating the relevant Green's 


functions and 


self-energies can 


be written in the following compact forms 






Gq'^ = G^ Q — TiQ i — So,/? , 


(7) 




G = Gj — Si — TiR , 


(8) 




/y—l _ 1 -y 


(9) 




" — "0 ^corr ■ 


(10) 



Here, the subscript i refers to an isolated molecule. Making use of the above equations, 
the Green's function matrix G of the molecule linked to electrodes needed to compute the 
current of Eq. ([T]) can be expressed in the following two equivalent forms 

G = Gj — So,L ~ So, /J — Sj, corr + ^Sj.corr ~ Scorr^ ; (H) 
G = Gj 0^ ~ Ti^corr — Tq^l — So,R — ^So,L " + So,/J — S/jj . (12) 

III. THE CONCEPT OF A SUFFICIENTLY LARGE EXTENDED MOLECULE 

The approach exposed below relies upon the concept of a sufficiently large extended 
molecule. We shall make the physical assumption that, even if, as a result of the molecule- 
electrode couplings, the electron correlations do not strictly remain confined within the 
molecule, they will be wiped out and eventually become negligible sufficiently far away 
from the molecule-electrode interfaces. Rephrasing, we shall consider physical situations 
where, although the intramolecular electron correlations can be strong, they rapidly decay 
beyond a characteristic length Cc and become negligible. It is the length £c that defines the 
size of the "sufficiently" large extended molecule. Being negligible beyond the sufficiently 
large extended molecule, correlations insignificantly affect the embedding self-energies which 
couple the extended molecule to electrodes (x = L,R), S^, — So,x ^ '^o,x- That is, the 
expression in the parenthesis of Eq. (|T2l) can be ignored if the extended molecule is sufficiently 



large. Importantly, this assumption is physically self-consistent: neglecting correlations 
outside the sufficiently large extended molecule implies that they are significant only inside it, 
"Scorr — '^i,corr, 1- G., the pareuthesis of Eq. (ITT]) also vanishes. The mathematical equivalence 
of these two physical conditions is guaranteed by Eqs. flTT]) and f[T^ . which show that the two 
parentheses vanish simultaneously, and yield a unique equation in the case of a sufficiently 
large extended molecule. Ultimately, the above physical assumption of a sufficiently large 
extended molecule is justified by just the problem that one wants to describe. If the opposite 
were true, the separation in uncorrelated left and right electrodes, and an extended (possibly 
strongly) correlated molecule would be impossible, and the problem would be no longer that 
of the transport through a single molecule. 

Practically, it is important that the extended molecule be not too large, such that accurate 
ab initio quantum-chemical calculations to correlated molecules of sizes Cc be feasible. As 
any approximation, the present approach also has a limited applicability. Importantly is 
however that — unlike other approximate schemes justified mathematically rather than 
physically — its limitation is physically transparent. The present approach cannot reliably 
tackle situations characterized by too large correlation lengths Cc- Fortunately, such cases 
are rare. We can only mention a sirigle case: the Kondo effect, extensively studied from 
various persepctives (e. g., Refs. I2CH24J and citations therein). As it is widely accepted, 
the Kondo anomaly in the zero-bias conductance is due to electron correlations related to 
coherent spin fluctuations within the Kondo cloud of linear size Ck ~ ^vp / [kBTx)- For 
typical Kondo temperatures Tk ~ ImK— IK and vp = 1.4 x lO^m/s (Fermi velocity in 
gold electrodes), Ck ~ 10^ — 10^ nm, a size which is well beyond the possibilities of any 
nontrivial ab initio study. As discussed recently [2^, Kondo-related properties obtained by 
treating exactly clusters of (much) smaller sizes, Cc <^ Ck, should be considered with very 
special care. 

We do not attempt to give a general estimate of the size Cc-, a hardly feasible task for real 
systems. Instead, we shall benchmark the approach based on a sufficiently large extende 
molecule for a non-trivial correlated model. On this basis, we argue that this approach can 
provide results, which are relevant for the transport in correlated molecular systems. 

So, if the extended molecule is sufficiently large, one can reliably determine the matrix 
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Green's function G of the embedded molecule as 



G-^ ^ GTl - 



G-^ ^ G7^ 



■'0,L 



■'0,R ! 



(13) 
(14) 



provided that one is able to accurately compute the matrix Green's function Gi of an iso- 
lated, but still sufficiently large extended molecule. In view of the physical motivation 
behind it, we shall refer to the approximation involved in Eqs. f lT3|) and f|T4|) as the approx- 
imation of a sufficiently large extended molecule. This concept plays a central role for all 
the subsequent considerations of the present paper. 

According to Eq. (IT^ . to compute the full matrix Green's function G of a sufficiently large 
extended molecule, in addition to Gi, one only needs Tx and Eq.x {x = L,R). More precisely, 
what one only needs is the retarded self-energy So,a;) because all the other quantities can 



be expressed in terms of it as follows [l6|, l25 | 



(15) 



-t (1 - /,) 



Consequently, the hard part to calculate the current is to obtain the matrix Green's 
function Gi of the isolated extended molecule. [As shown below, it even suffices to compute 
the retarded block GI, cf. Eqs. (117]) and (Il6]).] By inserting Eqs. ©, (E]), (ED, and (Il5]) 
into Eq. ([1]), one can express the steady-state current through a sufficiently large extended 
molecule as 

/ ^ ^ 1 deTr [hie) - fnie)] rL{e)G^{e)rR{e)G''{e) . (16) 

In the absence of correlations G^'"' — t- Gq'"', Eq. flTBl) . which reduces to the familiar Landauer 
brmula, is an exact result probably first deduced in Ref. and often discussed later, e. g. 



16 



18 



25 



27 



. Although similar to the uncorrelated case, Eq. ( [T6l) does account for 
electron correlations, because (G^'°')~^ — (Gq"') ~ —Hi^corr 7^ 0. It becomes similar 
to Eq. (10) of Ref. ilO, deduced for proportionate coupling, if one approximates S ~ So 



25 



if the 



there, in the spirit of the present paper. It is also similar to Eq. (17) of Ref. 
correction factor used there is set to unity (A = 1), that equation coincides with the present 
Eq. (|T6i) . The physical motivation of the coincidence of these two equations is clear: as seen 
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in the definition of A [cf. Eq. (13) of Ref. |25|, this factor accounts for the renormahzation 
of the electrode self-energies and should approach unity if the extended molecule becomes 
sufficiently large 29 |. 

Eq. fITB]) seems odd for the correlated it only contains the retarded Green's 

function [notice that = {G'^Y] but not the lesser Green's function [cf. Eq. ([T])]. However 
there is no contradiction. needed in Eq. ([1]), can be expressed in terms of G^ only if 
the bare embedding self-energies Ti^fi of Eq. (|T5l) are used instead of the dressed ones Ti^ in 
Eq. ([3]), and this is legitimate only for a sufficiently large extended molecule. For smaller 
sizes, both (nonequilibrium) G^ and G^ are needed to compute the current. 

The retarded Green's function block G^ of the sufficiently large extended molecule cou- 
pled to electrodes can be obtained from that of the isolated one G^ from a Dyson equation 
similar to Eq. ([2]) 

{G^-' ^ (G:)-' - SS,^ - (17) 

In view of Eq. (ITTI) . the retarded Green's function of an isolated but sufficiently large 
molecular cluster represents the crucial quantity, which is needed to compute the current 
from Eq. (fT6|) . For real molecules, G^ can be obtained, e. g., within diagrammatic schemes 
to approximate the self-energy S^^„^^. One of these schemes are the ADC- approximations 
widely employed for accurate electronic structure calculations Another scheme is 

based on the GW approximation, which was proposed long time ago for metallic solids [11 1 



and used in several recent studies of small molecules [33|, |3J| and nanotransport 5l-ll0|. It 
is certainly of interest to attempt to use such schemes for real systems, where exact results 
are not available, but applying them here in conjunction with Eqs. f|T3|) and f|T6|) is less 
useful conceptually. They represent certain mathematical approximations, which can be 
justified by comparing their predictions with experiments, but their validity in physically 
unclear. By constrast, our approximation possesses a precise physical content. Electrons can 
be strongly correlated in the molecule, but if interactions responsible for these correlations 
remain localized within a restricted spatial domain, the larger the extended molecule, the 
more accurate are the results; they should become exact in the limit of an infinite extended 



molecule 



25 



35|. 



Our present goal is to document the approximation of a sufficiently large extended 
molecule, and, to facilitate the interpretation of the results, we do not use any further 
approximation. Therefore, we shall employ a simple model of a correlated molecule, which 
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we solve exactly within a full CI numerical calculation. The exact retarded and advanced 
Green's functions of the isolated extended molecule will be computed with the aid of the 
Lehmann representation (case of zero temperature) 19 1 



r^wxi \r^a(.^* (0| aj |A) (A| 4 10) , (0| (y\ 4 |o) 



+ ^0^ 



e + e\r + «0 



Here, and are second-quantized operators that annihilate and create electrons at the 
sites ( "atomic orbitals" ) j and k within the extended molecule, whose ground state is denoted 
by |0), e^J" and stand for the v^^ ionization energy and A*^ electro-affinity, and {|z^)} 
and {|A)} represent the corresponding complete sets of eigenstates. 

IV. THE MODEL 



To be specific, we shall consider below a one- dimensional discrete model of a two-terminal 
setup known as the interacting resonant level model (IRLM) 3614391] . The physical system 
it describes is a molecule modeled by a single level of energy a attached to the left (L) and 
right [K) electrodes via nearest-neighbor electron hopping characterized by the resonance 
integrals /S^, and respectively. The semi-infinite electrodes are kept to fixed chemical 
potentials /iL,i?, and contain noninteracting electrons hopping between adjacent sites (reso- 
nance integrals C,l,r)- By ignoriiig the electron spin, subtle coherent fiuctuations (like those 
responsible for the Kondo effect 3, Isil) are disregarded, but the Coulomb interactions Jl,r 
at the molecule-electrode contacts, which are included, induce nontrivial electron correla- 
tions. They bring about an interesting physical behavior that attracted theorists' attention 
recently 36l-l39| . The important role played by interaction at the electrode- molecule contacts 
also made the object of recent experimental investigations 40 1. 

To avoid problems related to the Fermi energy alignment it is helpful to keep a constant 
band filling factor; it will be chosen 1/2, which amounts to consider a number of sites twice 
the number of electrons. We shall split the total Hamiltonian H of this system to explicitly 
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write the part Hm that specifies the extended molecule 



H 



Hl + Hr + Hlm + Hmr + H^ 



Hl = 

Hr = 

Hlm 
Hmr 



— oo — oo 



1=-Ml-1 
+00 



1=-Ml-1 
+00 



(19) 



and 



-Ml+1 



H 



M 



fiL Y Yl {^l^i-i + ^• 

l=-l 1=^1 

A/fl Mr-1 

r=l r=l 



(20) 



+Jl I a^ao 



In the above equation, the extended molecule comprises, besides the "molecule" itself, repre- 
sented by a single site (labeled by 0), parts of the left and right electrodes (sites —Ml, . . . , — 1 
and 1, . . . , Mr, respectively). Importantly, the energy of the molecular level a can be tuned 
by direct molecular orbital gating, as demonstrated by a remarkable very recent exper- 



iment 



4l| . or electrochemical gating 42|-|45| in a way that is simi 



"plunger" gate electrode in quantum-dot nanoelectronic devices 46|, |47|] . Unlike in quantum- 



ar to the usage of a 



dot based single-electron transistors, where the relevant a-range is of the order of a few meV 



471 ]. in the aforementioned experiment 



4l| the energy of the molecular level a could 



be varied, remarkably, within much broader ranges: 1.1 eV ^ £f ~ ^ 1.9 eV for single- 
electron transistors based on 1,8-octanedithiol and O.ieV ^ sp — C( < 1.8 eV for those 
based on 1,4-benzenedithiole by varying the gate potential Vg within —3.3 V < Vq < and 
—3 V < Vg < 3V, respectively. In the presently considered case of half filling, the number 
of sites of the extended molecule N = Ml + Mli + 1 must be even. This implies that Ml 
is odd and is even or vice versa, and therefore it is impossible to consider perfectly 
symmetrical extended molecules. The numerical results presented below correspond to the 



10 



most symmetrically possible situation Ml = Mr ± 1, but the conclusions of the present 
paper are not qualitatively affected by this choice. 

For one- dimensional models with nearest-neighbor hopping (the class to which the above 
model belongs), the contact "surfaces" between the extended molecule and the left and 
right electrodes reduce to two points, which are located at —M^ and +Mji, respectively. 
Therefore, the self-energies Sq and Sg^i? needed to calculate the current of Eq. (fT6|) possess 
a single nonvanishing elements, namely (Sj;^^)_^^^ = and (SJJ^jj)^^^^^^^ = S^j. 

For model ( !T9|) . ^^(e) = f^'^g^^ie) where is the retarded Green's function at the ends of 
semi-infinite electrodes coupled to the extended molecule. By using the expression of gl{e) 



2T\ . one straightforwardly obtains 



s.(5) = ^(5-/i,.)-^r.(^), (21) 



r.(e) = §V4e-(5-/iJ2 ^(2^^ - |5 - /X.I) , (22) 

where 9 is the step (Heaviside) function. The above expressions, along with the Dyson 
equation for the retarded Green's function deduced from Eq. IHM 

complete the formulae needed to compute the current (fT6ll . 



V. NUMERICAL RESULTS FOR THE LINEAR CONDUCTANCE 

We note that, although Eq. ( |T6l) as well as the subsequent equations for the nonequi- 
librium Green's functions are valid for an arbitrary voltage, in the presentation of the nu- 
merical results of this section we shall restrict ourselves to the case of linear conductance 
G = dl /dV\v=o- For completeness, we give its expression for the case of a sufficiently large 
extended molecule, which can be obtained by taking /i^ ^ = and V — )■ in Eq. (fT6|) 

G/Go ^ rL{eF)\GLM,,Mni^F)\'rn{eF) , (24) 

where Gq = e'^/h is the conductance quantum. For the sake of simplicity, we shall suppose 
that the electrodes are identical and therefore = = ^, Pl = Pr = /3, ^l{,£f) = 
^r{£f) = 2f3'^/C, = r, and = = J. Unlike in most of the studies on the IRLM or 
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similar models (see, e. g., Refs. Ua, I36M38I. l48l ). where ^ is chosen as the unit of energy, we 



shall use throughout /3 = 1. The reason is that we aim to compare the cases of organic 
and metallic electrodes among themselves, and their bandwidths (4^ within our model) 
significantly differ. On the other side, we do not expect substantial differences between 
the resonance integrals f3 corresponding to a molecule attached to organic and metallic 
electrodes, which should be of the order of a few electron volts. Basically, the same can also 
be expected for the Coulomb contact interaction J. 

Concerning the size N of the extended molecule we note the following. The main interest 
of the present paper is not the IRLM itself. Rather, we aim to unravel what one can 
learn from the results for this model for theoretical approaches to electron transport in 
real systems based on feasible ab initio quantum-chemical calculations able to accurately 
treat strong electron correlations. With the presently available computational resources it 
is unlikely that accurate calculations can include more than a few electrode layers. Within 
the present model, this would correspond to at most ~ 12. Therefore, it would make 
little sense to examine the results for the IRLM up to the largest size (A^ ~ 26), which 



can be tackled within full CI calculations by calculating the Green's 



unction of the isolated 



molecule from Eq. f|T8l) based on the Lanczos algorithm 2J, |48N53| , nor to carry out a 
finite scaling analysis of the numerical results to extract the limit N ^ oo. For the sizes of 
interest, the straightforward exact numerical diagonalization can be and has been employed 



to compute the exact Green's function Eq. ( jlSll . within full CI calculations. 
A. Case of organic electrodes 

To model a (chemisorbed) molecule linked to organic electrodes, one should consider 
that the resonance integrals ^ and /3 are of comparable magnitudes, and this results in large 
values of the hybridization parameter F = 2/3^/^. For illustration, in this subsection we 
shall examine the case where they are equal, ^ = /3(= 1). 

First, we want to benchmark the approximation of the sufficiently large extended molecule 
presented in Sect. [TTl against a well established exact result, namely the case where the 
molecular level is at resonance with the electrode Fermi energy {a = ep)- There, similar to 
ihe cases of the noninteracting resonant level (J = 0) and of the single-electron transistor 



20|, |2l|, |47| , the conductance of the IRLM reaches the unitary limit G = Gq. Physically, this 
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results from the fact that (irrespective of the values of ^, /3, and J) for = a = hr all sites 
are half occupied, rij = (0|ajaj|0) = 1/2, and this yields, via Friedel's rule G = Gosin^(7mo) 
a perfect transmission {G = Gq). Numerically, the result G = Gq on resonance was 



obtained by means of time- dependent density matrix renormalization group calculations 



37|. 



Our results for the on-resonance case at /3 = ^ are presented in Fig. [2^. As visible 
there, even the smallest possible extended molecule [N = 4) represents a "sufficiently large" 
extended molecule, enabling to accurately reproduce the exact result. Indeed, even for the 
strongest coupling shown in Fig. [2^ (J = 10, which is in fact unrealistically strong), the 
numerical results deviate from the exact value G = Gq by ~ 13%. As expected from the 
physical analysis backing the approximation of a sufficiently large extended molecule. Fig. [2] 
reveals that the deviation from the exact result diminishes with increasing sizes; for the next 
smallest size (A^ = 6), the error is at most 5%. 







N=4 




N=6 




N=8 




N=10 




— N=12 




N=I4 






FIG. 2: Conductance G computed numerically from Eq. ()24p for several sizes in units of quantum 
conductance Gq = e^//i plotted versus: (a) Coulomb contact repulsion J at resonance (a = £f), 
and (b) molecular level energy a for ^ = 1, J = 4. The sizes of the extended molecules are given 
in the legend. The thin dashed line in panel (b) corresponds to a vanishing Coulomb interaction 
(J = 0). 



Because G = Gq holds on resonance even for the uncorrelated case (J = 0), where the 
results are anyway independent on the size of the extended molecule (cf. Sect. I), one might 
suspect that the very weak size-independence displayed in Fig. [2^ does not demonstrate the 
validity of the approximation of a sufficiently large extended molecule, but rather that the 
electron correlations themselves are weak. To illustrate that this is not the case, we present 
in Fig. [2b results for the conductance out of resonance {a ^ Sp) and realistic, moderately 
strong Coulomb strength J = 2r = 4 (and equal to the electrode bandwidth 4^ = 4). 



13 



They confirm the above finding, namely that the size-independence is indeed very weak. In 
addition, the conductance for J = 4 can be compared with that for the noninteracting level 
(J = 0) (thin dashed line in Fig. the substantial difference between them indicates a 
significant electron correlation effect on the electric transport. To exemplify, we note that 
this Coulomb interaction diminishes the half-width of the curve G = G{a) roughly by a 
factor two (from F = 2 to 0.8). AtQ; = £:i7'±F(F = 2/3^/^ = 2), i. e. at the half maximum 
in the uncorrelated case, the electron correlations suppress the conductance by a factor ~ 5. 
Importantly, for situations even farther away from resonance (and this is the usual case. 



4ll ) , Fig. shows that the suppression 



including also that of the recent experiments of Ref. 
factor becomes considerably larger. 

On the basis of the results presented in this subsection, we claim that the correlated 
electron transport through a molecule attached to organic electrodes can be accurately 
obtained by carrying out microscopic calculations for small extended molecules; besides the 
(true) molecule of interest, they should only include a reduced number of electrode layers. 



B. Case of metallic electrodes 



In comparison with organic electrodes, metallic electrodes are characterized by a larger 
bandwidth (4^). To model (physisorbed) molecules weakly coupled to metallic electrodes 
within the framework of the present paper, one should consider values of ^ larger than in 
Sect. IVAl 

In Fig. [3l we present numerical results on resonance for two values of the molecule- 
electrode resonance integral, ^ = 10 (F = 0.2) and ^ = 5 (F = 0.4), and several sizes 
of the extended molecule. The trend that, as the extended molecule becomes larger, the 
linear conductance approaches the unitary limit {G — Go) can be seen in both panels of 
Fig. [3l So, the basic physical argument behind the approximation of a sufficiently large 
extended molecule is also confirmed by these numerical results. However, unlike in Fig. [2^, 
the deviations of the numerical values of G from the exact value Gq visible in Fig. [3^ and b 
are significant and deserve further analysis. 

Figs. IHK and b show that the departure from the unitary limit is larger for a stronger 
Coulomb contact strength J. This behavior can be easily understood. The stronger the 
electron-electron interaction, the more important is the renormalization of the contact self- 



14 



energies "Ex from the bare self-energies I]o,£c- This renders more significant the quantity 
in the parenthesis of the r.h.s. of Eq. f lT2|) . which was neglected to derive the approximate 
Eqs. f ll6p and fl24p . To make those parentheses negligible, one has to reduce the renormal- 
ization at the ends of the extended molecule due to electron correlations, and for this one 
has to sufficiently increase its size A^. 

As illustrated by Figs. [3^ and b, for larger ^'s, the unitary limit can be reasonably 
reproduced by using the smallest possible size = 4 only for weak Coulomb interactions 
of the order of the bare hybridization J ~ F. Unfortunately, such Coulomb strengths are of 
little interest. They are too weak, and the correlation effects they bring about are negligible. 
The counterpart of Fig. the linear conductance G = G{a) for J ~ F is not significantly 
different from that for J = 0. For ,^ = 10, a Coulomb strength J = 2F = 0.4 solely yields 
a slight change of the half-width at half maximum from F = 0.2 to 0.22, and for = 5 and 
J = 2F = 0.8 the change in the half- width is only a bit larger (from F = 0.4 to 0.485). 

The impact of the stronger but still realistic Coulomb contact interaction of Fig. |2b 
(J = 4) is also significant at larger ^. For ,^ = 10 it increases the half-width at half maximum 
of the curve G = G{a) from F = 0.2 to 0.43, whereas for ^ = 5 the increase is from F = 0.4 to 
0.80. However, the corresponding values of J/F = 20 and J/F = 10, respectively, are much 
larger than that of Fig. |2]d (J/F = 2). At these values of J/F, the curves of Figs. [3^ and b 
exhibit a notable A^-dependence. However, for the smallest extended molecule (A^ = 4) the 
computed conductance still amounts a significant fraction of the exact conductance (~ 45% 
and ~ 55%, respectively). For A^ = 8, which corresponds to include three and four layers in 
the left and right electrodes, these fractions are increased to ~ 67% and ~ 80%, respectively. 
In view of the state-of-the art in the field of molecular electronic transport, characterized by 
large differences between theory and experiments, these values suggest that, although less 
accurate than for organic electrodes, achieving a semi-quantitative description based on ab 
initio quantum-chemical calculations is possible even in the case of metallic electrodes with 
moderately strong Coulomb contact repulsion. 

VI. A SIMPLE RENORMALIZATION SCHEME 

As alternative to the approximation of a sufficiently large extended molecule, one can 
start from the exact Dyson equation ([8]) and attempt to use renormalized electrode self- 
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FIG. 3: Conductance G in units of quantum conductance Go = e^/h plotted versus contact 
interaction J at resonance a = Si? for ^ = 10 and ^ = 5 computed using extended molecules with 
N = 4,6, 8, 10, 12, 14, 16 sites (values increase upwards). 

energies Sa, instead of So,a;; by devising certain diagrammatic approximations. We shall 
not pursue this line here, because as noted in Sect. [TTl the validity of such approximations 
is usually not transparent physically. We present instead a simple renormalization scheme, 
which turns out to work surprisingly well for the IRLM. 



Inspired by the rather general framework of the Landau Fermi liquid theory [5J], we shall 
simply suppose that the aforementioned renormalization can be accounted for by means of 
a multiplicative e-independent factor F, G'^{e)\exact — FG'^{e). Then, the problem is to 
deduce this factor. For the IRLM at resonance, this can be done, because, as already noted 



in Sect. IV At the exact conductance is known, G\a=£p = Gq. This condition determines the 
needed factor at resonance for given values of A^, C, and J (remember that we keep /3 = 1), 
but F could dependent on a. Because G = G{a) reaches its maximum at a = ep one can 
admit that, at least not too far away from resonance, F only slightly depends on a and 
neglect this dependence altogether. 

We have computed and examined the curves of G = G{a) obtained within this renor- 
malization procedure. The results are very encouraging. The renormalized curves are much 
less size-dependent than the unrenormalized ones. The smaller the values of ^ and J, the 
weaker is the dependence on A^. For illustration, we present in Fig. H] for comparison the 
renormalized curves along with those unrenormalized in the rather unrealistic situation of 
very large values of J = 10 and ^ = 10, which corresponds to a very large ratio J/T = 50. 
The fact that, with increasing A^, the renormalized curves tend to become size- independent 
much faster than the unrenormalized curves is clearly visible even in this less favorable case. 
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FIG. 4: Normalized conductance G/Gq and several sizes A^. The lower group of curves (thin lines, 
small points) is computed by means of Eq. ()24p . and the upper group (thick lines, big points) are 
obtained within the renormalization scheme described in the text. Notice that, unlike the strongly 
size-dependent unrenormalized curves, the renormalized curves are weakly size-independent even 
for the rather large values J = 10 and ^ = 10 used, which are less favorable for the renormalized 
procedure employed. 

So far, we have considered the dependence on the location of molecular level a relative 
to the electrode Fermi energy. We also employed this renormalization scheme to study the 
J- and ^-dependence of the linear conductance away from resonance a ^ Sp. (Remember 
that G\a=sp = Go for all ^'s and J's.) Results of this kind are depicted in Fig. [5l which 
includes both unrenormalized and normalized curves (thin and thick lines, respectively). 
They confirm that the size-dependence is stronger at larger ,^'s, as already mentioned above. 
In addition, they reveal the interesting fact that the dependence of the Coulomb contact 
interaction J can be non-monotonic, as is the case of other previously published results 



37( 1 . As a general trend, we found that by increasing J the conductance G{J) gradually 
increases from the value G\j=o, reaching its maximum at a Coulomb strength Jc ~ ^, and 
decreases afterwards. This trend can be seen by inspecting Fig. [5l which also shows that the 
wider the electrode bandwidth (4^), the more pronounced is the conductance enhancement 
G(J,)/G(J = 0). 



VII. DISCUSSION 



Within the NEGF-DFT approaches, the impact of electron correlations on the Green's 
function of the isolated extended molecule is accounted for approximately, Gi\DFT ~ Gi, by 
employing certain exchange-correlation (xc-)functionals. To determine the Green's function 
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FIG. 5: Normalized conductance G/Gq plotted versus J/T for ^ = 10,5,3.33,2. All the curves 
were computed for the out-of-resonance situation \a — £f\ = T/2, which corresponds to G/Gq = 
0.8 at J = 0. Notice that, the smaller the ^-value, the stronger the size-dependence of thin 
curves computed from Eq. (|24|) . and the larger the departure from the practically A^-independent 
renormalized curves, depicted by thin lines. 



G\ DFT = ( G^i |dft ~ So,i — So,i? ) of the extended molecule coupled to electrodes, 
those approaches employ the bare self-energies So,i and So, it- Because the isolated extended 
molecule is modeled by an effective Kohn-Sham effective Hamiltonian, the underlying picture 
is a one-particle picture, and therefore the Landauer-type formulae (fT6|) and fl2^ with & 
replaced by G'^\dft follow as "exact" results (i. e., without any other approximation as the 
DFT-approximation itself). 

We emphasize that Landauer-type expressions like those Eqs. ( 1T6|) and ( 12^ are not valid 
in general for correlated electron systems. As explained in Sect. [Tll they represent justi- 
fied approximations in the correlated case only when sufficiently large extended molecules 
are considered. The advantage of the approach based on this concept is its clear physical 
content. Eqs. (fT6|) and (!24|) hold for electron correlations that can be strong within the 
molecule but are negligible at the ends of a sufficiently large extended molecule. This is the 
reason we employed in the Dyson equation the bare self-energies So,l and So,it, as in 
the NEGF-DFT approach. However, at variance with the latter, our approach accurately 
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accounts for molecular electron correlations, which are included via the Green's function 
of a sufficiently large extended molecule. Because our method is based on the exact Green's 
function (whose counterpart for real molecules is a Green's function computed by accu- 
rately including electron correlations within appropriate ab initio methods), it is clear that 
it will provide more reliable results than NEGF-DFT's. 

To give support to the intuitive physical reason underlying the concept of a sufficiently 
large extended molecule, we have presented numerical results deduced within a non-trivial 
correlated model. Based on them, we argue that the present approach is able to provide 
results, which are interesting for the transport in correlated molecular systems. To this 
aim, the essential prerequisite is the ability to accurately compute the Green's function 
of the isolated extended molecule. Fortunately, this can be regarded as a feasible com- 
putational task. Such Green's functions were already obtained by means of the Lanczos 



algorithm within the ADC 30|-|32|. The information contained in G[ [ionization energies 



and electro-affinities, and pole strengths, see Eq. ( ITSl) ] can be directly compared with ex- 
perimental results, which can be obtained from appropriate experimental investigations of 
the isolated molecule. Encouragingly, accurate ADC schemes were successfully applied for 
correlated medium size molecules (e. g., benzene), comparable to those of interest for molec- 
ular electronics |[4£|. Importantly, although these methods are much more accurate than 
the currently considered GW, the computational effort is not much larger; it scales as Af^ 
within ADC(2) and as ~ JV^ within GW, where AC is the number of atomic orbitals 55 |. 
The difference between the present approach and the usual NEGF-approaches to cor- 



related transport [SHlOl. l28| can be discussed by inspecting the parenthesis of the r.h.s. of 
Eq. ( fTTj) . In our case, this parenthesis becomes negligibly, and the difficulty to acount for 
correlations is to accurately tackle a sufficiently large extended molecule. If we considered 
the true (small) molecule and neglected that parenthesis, we would make an approxima- 
tion similar to the so- 
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Then, the current 



would be (rather poorly) approximated by Eq. ( IT6ll |28l |. This amounts to take = 3 
59[, which as seen by inspecting Figs. El HI and [5l is not a good approximation. Our ap- 
proach goes beyond the elastic cotunneling approximation, it accounts for all higher-order 
processes within the extended molecule. Other NEGF-approaches to correlated transport 
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(which are also beyond the elastic cotunneling) adopt a different standpoint: 
they consider the true (small) molecule linked to electrodes and the evaluate the aforemen- 
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ioned parenthesis (which does not vanish) by certain Keldys 



19|. Recent studies often employ the GW approximation 



1 diagrammatic approximations 
SMlO]. Results obtained within 



the GW approximation are certainly of interest [5|-l8|, although limitations and artefacts are 
also significant B Q. MethodOcca,,, the weaU pent of U„s app.oacK is the fact that 

its physical content is hard to understand. Thus, it is difficult to assess a priori its validity 
for a real case. The attempt to benchmark the GW approximation against exact results 

nn 

for transport led to ambiguous conclusions [9|, |10[ ■ Within this framework, both important 
effects (i. e., intramolecular correlations and molecule-electrode coupling) are treated in an 
approximate manner. Consequently, it is difficult to conclude which of them is better or 
worse approximated. It makes more sense to benchmark the GW approximation for those 
isolated extended molecules, wherein electron correlations play an important role, and which 
are also employed in transport experiments. Recent studies in this direction are significant. 



but only considered atoms and very small molecules so far [33|, |34 1 . 

In Introduction, we mentioned that the results for transport should be partition inde- 
pendent. Obviously, within the present approach, where correlations in distant parts of 
electrodes are neglected, the partition independence means and is demonstrated by the fact 
that, by making the extended molecule larger than a certain size Cc, the results do not 
(notable) change. Regarding the size- (in) dependence of the results of other approaches to 
correlated molecular transport, attention should be paid to the following important aspect. 
The fact that calculations accounting approximately for electron correlations (like DFT's or 
GW's) yield results that saturate with increasing N, does not necessarily imply that they 
are physically meaningful. The saturation could be a mathematical artefact, because such 
approximative methods become poorer at larger sizes and cannot reliably describe subtle 
correlations basically localized within the molecule in a system including larger and larger 
parts of uncorrelated electrodes. Correlation effects at larger sizes are likely underestimated, 
and this renders a correlated system resembling more an uncorrelated one, wherein the size- 
dependence is absent (cf. Sect. [11), and hence the aforementioned saturation. On the contrary, 
if molecular correlations can be accounted for by means of an accurate Green's function G\ of 
an isolated and sufficiently large extended molecule, size-independent results can be trusted: 
they indicate that correlations are negligible at the interfaces between electrodes and the 
extended molecule. 

The present numerical results indicate that the extended molecules to be employed for 
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carrying out accurate ab initio quantum chemical calculations need not be very large. Partic- 
ularly attractive is the case of organic conducting electrodes, for which only a few electrode 
layers have to be included (cf. Sect. IV Ap . Another advantage is that, without metal atoms, 
ab initio calculations are less demanding. For metallic electrodes, more electrode layers 
need be included (cf. Sect. IV Bp . but one can still hope to a semi- quantitative description. 
Alternatively, simpler (tight-binding) models can be utilized to include metallic electrode 
layers and their coupling to the molecule as well, similar to the first four terms of Eq. fl20p . 
Regarding the case of metallic electrodes, it is still worth emphasizing the following aspect. 
Experimental data indicate that electric conduction through molecules typically occurs by 
non- resonant tunneling, characterized by a rather lar ge e nergy offset {ep — EhomoI (p-type 



conduction) or Ielumo — £f\ (n-type conduction) [4l|, |60|, which represents the counterpart 
of the quantity |a — ei?! in the presently employed model. This fact is important because, 
based on our calculations, we can state that the size dependence of the transport properties 
is particularly pronounced near resonance (|a — e^l ~ 0), but becomes much less significant 
away from resonance. This feature, which is visible if one inspects the unrenormalized curves 
of Fig. m is very encouraging: one can hope to obtain quantitatively correct results by using 
extended molecules with reasonable sizes even in the case of metallic electrodes. 

There is a general claim (e. g., Ref. l6ll) that whether the electron transport is correlated or 
not solely depends on value of the ratio between the characteristic Coulomb electron-electron 
interaction and the hybridization (in our case J/T). By comparing the results of Sects. IV Al 
and IV Bl among themselves, we found limitations of this rule of thumb 611] • Based on these 
results, one can conclude that the ratio J/T alone is not sufficient to distinguish between 
significant and insignificant correlation effects on the transport. For large ^'s, significant 
effects arise at Coulomb strengths J substantially larger than F, while for C, ^ P values 
J ~ F suffice. 

As mentioned from the very beginning, the main interest of the present work was not to 
analyze the physical properties of a simple model (i. e., IRLM), to which all our numerical 
results refer. Rather, we attempted to extract from them information, which is useful from 
the perspective of real systems of interest for molecular transport. Nevertheless, in spite 
of its simplicity, we emphasize that the IRLM incorporates a significant physical aspect, 
namely the Coulomb repulsion at the contacts. 

Concerning the physical relevance of the employed IRLM, we also want to draw the 
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attention on the following fact. To give support to the robustness of their theoretical results, 
many investigators mentioned that reasonable changes of the molecule-electrode distance d 
do not notably alter the calculated transport properties. Such a behavior can indeed be an 
indication for robust results, but it could also have a different origin. Within the IRLM, an 
increasing d makes both parameters /3 and J to decrease. The decrease of /3 always causes a 
conductance reduction. As already noted and visible by inspecting the renormalized curves 
of Fig. [5l depending on the range, a decrease in J could either suppress or enhance the 
conductance. Noteworthy, this effect brought about by a modification in J can be observed 
on the renormalized curves. For realistic values of J(~ 1), the unrenormalized conductance 
increases with decreasing J. So, rather than indicating robust results, the aforementioned 
behavior could be related to the fact that the size of the extended molecule employed in the 
calculations is too small, and for that size there is a compensation of two opposite effects 
induced by the /3- and J-variations. 

Of course, we do not expect that in all cases the largest extended molecule that can 
be tackled is sufficiently large to justify neglecting the influence of electron correlations 
on the self-energies at the contacts. In Sect. IVIl we considered a renormalization scheme 
appropriate for the IRLM, for which the exact conductance on resonance is known. To 
avoid misunderstandings, we emphysize that this scheme represents an extra renormaliza- 
tion, which solely accounts for the fact that the parentheses of Eqs. ffTTl) and f|T^ do not 
exactly vanish. The important renormalization (the exact intramolecular correlations of the 
isolated extended molecule) is already incorporrated in G[. Simple schemes similar to our, 
which underline the renormalization considered in the Landau Fermi liquid theory, work 
pretty well even in Kondo context (see, e. g., the supplementary material of Ref. |62| ). Cer- 
tainly, the renormalization scheme of Sect. |Vl] cannot be applied in an identical form in 
general. However, we think that the possibility to renormalize the Green's function G'^{e) 
of the embedded molecule by an e-independent factor, which furthermore turned out to 
be (practically) independent on the model parameters, is interesting and deserves further 
attention. This (almost) constant factor can no more be determined from the unitary limit, 
as done in Sect. IVIl but one may think to determine it by usin g ge neral properties of the 



Green's functions (e. g., sum rules or asymptotic behavior [19|, l62|). The speculation on 
the existence of a (nearly) constant renormalization factor is also appealing from a different 
perspective. In the cases where the renormalization effect can be accounted for by an almost 
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constant multiplicative factor of the Green's function, the renormalized and unrenormalized 
curves for current or linear conductance roughly differ by a nearly constant scaling factor 
[cf. Eqs. ( ITB]) and (^^]- It is tempting to relate this to the fact that, in spite of substan- 
tial quantitative discrepancies, this is just the basic difference between the calculated and 
measure curves observed in some cases. 

VIII. CONCLUSION 

In the present paper, we have presented a theoretical approach to correlated molecular 
transport relying upon a clear physical concept, namely that of a sufficiently large extended 
molecule. The analysis of the present numerical results obtained for a simple but non-trivial 
correlated case has indicated that this approach can be legitimately applied to the molecular 
transport in real systems. Computing the Green's functions within accurate quantum chem- 
ical calculations appears to be a feasible task for the isolated medium-size molecules needed. 
This is particularly true for narrower-band (organic) electrodes, but wider-band (metallic) 
electrodes can still be tackled especially in the most frequent situations of non-resonant 
tunneling. To remedy the stronger size-dependence for wider-band electrodes closer to res- 
onance, we have suggested that simple renormalization procedures inspired by the Landau 
Fermi liquid theory can work, because it is unlikely that electron correlations in the systems 
of interest for molecular transport could invalidate the latter. Although our main interest 
was not the IRLM, we still note that the counter-intuitive theoretical prediction, that the 
contact electron-electron interaction could enhance the conductance, is interesting and can 
make the object of a parallel experimental investigation. 
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